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ABSTRACT 

Triplex  arithmetic  is  a variant  of  interval  arithmetic  in  which  a "main" 
value  is  computed  in  addition  to  the  endpoints  of  the  containing  interval. 
The  "main"  value  is,  in  general,  the  value  that  would  have  been  comDuted  had 
properly-rounded  real  arithmetic  been  used  to  compute  the  results  of  the  cal- 
culation; it  may,  in  some  sense,  be  regarded  as  the  "most  probable"  value  of 
the  result  of  the  calculation*  The  endpoints  of  the  interval  define  the 
worst-case  range  of  the  possible  values  of  the  computation* 

This  report  describes  a FORTRAN  implementation  of  triplex  arithmetic  in 
both  single  and  multiple  precision*  The  package  described  in  this  report  is 
designed  to  be  used  with  the  AUGMENT  precompiler;  this  makes  triDlex  arithme- 
tic particularly  easy  to  use* 
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SIGNIFICANCE  AND  EXPLANATION 


In  some  computations,  it  is  essential  to  determine  rigorous  bounds  on  the 
computed  results*  Interval  arithmetic  can  be  used  for  this  purpose;  the 
endpoints  of  the  computed  intervals  Rive  rigorous  bounds  on  the  range  of  the 
result  values.  One  of  the  difficulties  with  interval  arithmetic  is  that  the 
intervals  may  become  so  wide  as  to  be  meaningless  in  some  problems*  One  might 
elect  to  take  the  midpoint  of  the  result  interval  as  being  an  approximation  to 
the  value  of  the  result,  but  this  may  not  be  entirely  satisfactory,  either  — 
for  example,  if  the  result  interval  were  [-infinity,  infinity]  the  midpoint 
would  be  computed  to  be  zero  — but  this  is  hardly  a meaningful  result,  since 
the  midpoint  could  be  any  real  number. 

In  triplex  arithmetic,  a "main"  value  is  computed;  this  value  is  the  re- 
sult that  would  have  been  obtained  had  properly-rounded  real  arithmetic  been 
used  in  the  computation*  This  gives,  in  some  sense,  a "most  probable"  value 
for  the  result* 

In  this  paper,  the  authors  describe  a program  package  which  performs  cal- 
culations in  triplex  arithmetic.  The  package  is  written,  insofar  as  practica- 
ble, in  ANSI  standard  FORTRAN,  and  it  is  designed  to  be  used  in  conjunction 
with  the  AUGMENT  precompiler  for  FORTRAN  to  afford  the  user  a relatively 
painless  way  of  carrying  out  computations  in  Triplex  arithmetic* 


The  responsibility  for  the  wording  and  views  expressed  in  this  descriptive 
summary  lies  with  MRC,  and  not  with  the  authors  of  this  report* 


TRIPLEX  ARITHMETIC  FOR  FORTRAN 


Klaus  Boehmer  and  J.  M.  Yohe 


Interval  arithmetic  is  a means  of  obtaining  rigorous  bounds  on  the  re- 
sults of  computations  by  carrying  both  upper  and  lower  bounds  through  the  com- 
putation. This  method  was  first  formalized  by  R.  E.  Moore  [11],  and  is  dis- 
cussed in  detail  there  and  elsewhere  ([1],  [12]).  Briefly,  however,  the 
bounds  obtained  by  interval  computations  represent  worst-case  bounds  for  all 
facets  of  the  computation:  error  of  approximation  in  the  input  data,  error  in 
obtaining  the  input  data  (via  measurement,  for  example),  and  errors  introduced 
in  the  computation. 

One  of  the  difficulties  with  interval  arithmetic  is  that  the  intervals 
may  become  so  wide  as  to  be  meaningless  in  some  problems.  One  might  elect  to 
take  the  midpoint  of  the  result  interval  as  being  an  approximation  to  the  val- 
ue of  the  result,  but  this  may  not  be  entirely  satisfactory.  Let  us  suppose, 
for  example,  that  we  are  working  with  a 6-place  decimal  computer.  Using  in- 
terval representation,  we  would  write 

JI=  C3. 14159,  3.14160], 

Note  that  11=  3*1415926535*. • belongs  to  the  interval  JI • II  is  a pretty  good 
approximation;  however,  we  would  find  that 

n2  = [9.86958,  9.86966] 

--  a spread  of  8 in  the  last  position.  If  we  taka  the  midpoint  of  this,  we 
have 

tt2  = 9.86962 

— which  does  not  compare  favorably  with  the  correctly  rounded  value  of  it2  = 
9.86960^....  We  would  in  fact  obtain  9.86959  if  we  were  to  apply  correct 
rounding  to  the  product  3*14159  x 3.14159*  The  reason  for  this  is  that  often 
the  round-off  errors  in  computations  cancel  one  another  in  a very  advantageous 
way,  so  a result  obtained  by  such  a correct  round-off  procedure  is  often  much 
better  than  the  endpoints,  or  even  the  midpoint,  of  the  computed  interval. 

Therefore,  it  has  sometimes  been  proposed  that  an  interval  be  reoresented 
as  a midpoint  together  with  an  error  term  e (this  is  known  as  the  centered 
form).  In  theory,  this  is  Just  as  good  as  the  representation  in  terms  of 
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endpoints;  in  practice,  it  is  not#  If  we  were  to  represent  n in  centered 
form,  we  would  need  to  write 

H = 3.14159  + .00001 

This  would  give  us 

n2  = 9.86959  ± .00007, 

generating  an  upper  bound  which  is  no  better  and  a lower  bound  which  is  worse. 
(We  acknowledge  that,  in  the  case  of  high  precision  calculations,  the  addi- 
tional error  introduced  by  this  method  is  not  so  serious,  and  considerable 
storage  economy  may  be  realized  by  using  the  centered  form.  This  was  dis- 
cussed by  Aberth  in  [1]  and  in  [2].) 

We  also  acknowledge  that  it  is  possible  to  gain  extra  precision  in  the 
Centered  form  by  eliminating  leading  zeros  in  the  error  term,  thus,  for  exam- 
ple, representing  II  by 

n = 3.14159  ± 2.65359  x 10-6. 

In  this  case,  we  would  obtain 

H2  = 9.86959  + 1.89449  x 10-5 

for  a low-order  spread  of  about  3.8.  This  is  more  accurate  than  standard  in- 
terval arithmetic,  but  requires  additional  computing  effort  which  might  not  be 
realizable  or  practical  on  a given  computer. 

Alternatively,  we  may  carry  three  numbers  through  the  computation:  UDDer 
and  lower  bounds,  as  for  intervals,  and  a main  value,  which  represents  the 
correctly  rounded  result  of  the  computation.  Using  this  scheme,  we  would  have 
the  triplex  numbers 

n*  = [3.14159,  3.14159,  3.14160], 


and 


n*2  = [9.86958,  9.86959,  9.86966], 

In  this  case,  the  main  value  is  very  close  to  the  correctly  rounded  value  of 
the  result.  Obviously,  it  need  not  be  so,  and  the  only  interpretation  that 
can  be  given  to  the  "main"  value  is  that  it  is  "probably"  "close"  to  the  cor- 
rect answer. 

This,  then,  is  the  motivation  for  triplex  arithmetic.  If  we  were  not 
concerned  with  finite  precision  calculation,  we  would  not  need  to  distinguish 
between  triplex,  the  centered  form,  and  interval  endpoint  representation  of 
data,  since  they  would  all  give  the  same  result.  In  the  context  of  fixed  fi- 
nite precision  representation,  however,  triplex  can  provide  us  with  additional 
insight. 
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The  theory  of  triplex  arithmetic  is  developed  extensively  in  [11]  and 
[12],  for  example,  and  we  refer  the  reader  to  these  sources  for  a more  compre- 
hensive discussion* 

The  first  complete  system  for  performing  triplex  computations  that  is 
known  to  us  was  the  Triplex-ALG0L-60  compiler  developed  by  Wipperman  and  oth- 
ers at  the  University  of  Karlsruhe  [15]  in  the  late  1960's*  A previous 
triplex  arithmetic  package  for  FORTRAN  was  developed  at  the  Mathematics  Re- 
search Center,  University  of  Wisconsin  - Madison,  by  Boehmer  and  Jackson  [2] 
in  about  192 * * * * 75;  this  package  comprised  a complete  single  precision  version 
(based  on  an  earlier  version  of  the  interval  arithmetic  package;  see  [4])  and 
a multiple  precision  version  with  limited  capabilities*  Both  versions  were, 
however,  heavily  dependent  on  the  UNIVAC  1100  architecture,  and  conseauentlv 
lack  the  portability  and  flexibility  of  the  present  package* 

In  this  paper,  we  describe  a complete  triplex  arithmetic  package  for  FOR- 
TRAN, based  on  the  interval  arithmetic  package  developed  by  the  first  author 
([18],  [I7])*  This  triplex  package  comes  in  two  versions:  a single  precision 
package  which  is  reasonably  portable  (see  [18]),  and  a completely  portable 
multiple  precision  version  based  on  the  multiple  precision  arithmetic  package 
of  Brent  [3]* 

Since  there  are  strong  similarities  between  this  triplex  Dackage  and  the 
interval  arithmetic  package  described  in  [I7],  we  will  not  reformulate  all  of 
the  presentation  given  there,  but  will  refer  to  that  paper  where  it  seems  ap- 
propriate* 

The  current  package  is  designed  to  be  used  with  the  AUGMENT  precompiler 
for  FORTRAN  (see  [5],  [6],  [7]),  and  indeed  is  itself  written  with  the  aid  of 
the  AUGMENT  precompiler*  This  is  discussed  more  fully  elsewhere  ([8],  [I7]), 
and  will  not  be  iterated  here.  Of  course,  for  the  multiple  precision  version, 
an  AUGMENT  interface  for  Brent's  package  must  be  used;  this  is  described  in 
[4]. 


2.  A formal  definition  of  triplex  rjurgfeecs.: 

A triplex  number  is  an  ordered  triple 

T :=  (tL,  tM,  tR)  with  tL  i tM  i tR  (2.1) 

We  call  tL  the  left  endec/ . tM  the  main  value,  and  tR  the  righl  entiPoj-lU  of 
the  triplex  number* 

We  denote  the  interval  defined  by  T to  be 

IT  :=  [tL,  tR]  (2.2) 

so  all  three  components  of  the  triplex  number  are  elements  of  the  interval  IT. 
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t 


If  tL  = tR  (=  tM  by  definition)  we  say  that  T is  degenerate*  Degenerate 
triplex  numbers  may  be  identified  with  the  corresponding  real  numbers;  hence 
the  reals  may  be  regarded  as  a subset  of  the  triplex  numbers* 

We  want  to  point  out,  however,  that  the  conversion  of  a real  number  £ on 
computer  in-  or  output  does  give  an  approximation  for  an  unknown  real  number 
t,  usually  without  any  information  about  how  well  or  how  badly  £ approximates 
t c R.  If  we  find  on  output  a degenerate  triplex  [£,  £,  £] , approximating  t, 
we  know  that  £ = t (see  below)*  So,  at  least  in  this  connection,  there  is 
much  more  information  in  a degenerate  triplex  number  than  in  the  corresponding 
real* 


In  finite  precision  calculations  the  situation  is  complicated  by  the 
need  for  rounding*  We  first  have  to  presume  that  rounding  is  performed  ac- 
cording to  the  mathematical  foundation  laid  down  by  Kulisch  [9]*  If  M is  the 
set  of  machine  numbers,  ri:R  ->  M is  said  to  be  an  optimal  rounding  if  and  onlv 
if 


x£  y =>  tf(x)  £ 4>(y) , and 


(2.3) 


d(x)  = x for  all  x e H 

Equation  (2*3)  implies  that  any  real  number  which  is  not  a machine  number  must 
be  rounded  to  one  of  the  two  machine  numbers  which  bracket  it: 

Lemma:  If  4 is  an  optimal  rounding,  and  if 

m^  :=  max  li  I I e M,  x < x)  < x < m2  :=  min  {x.  i £ e M,  £ >.  x} , 

then  d(x)  = m^  or  ri(x)  = m2  • 

Proof:  We  have  d(tn^)  £ d(x)  £ t (m2)  by  (2.3a)*  But  by  (2*3b),  d(m,)  = m^, 

and  ri(m2)  = m2,  so  m-^  £ <tf(x)  £ m2*  Since  ri(x)  e M and  M n [m^,  m2]  = 

Im-p  m2),  we  have  d(x)  = m^,  or  ri(x)  = m2*  [] 

Now  we  define  a correct  rounding  6 of  a triplex  number  into  the  set  TM  of  ma- 
chine triplex  numbers: 

V(x)  is  the  "optimal  downward  directed  rounding"'  i.e* 


( 2* '' ) 


V(x):=  max  {£.  ! i e M,  x £ x)  for  all  x e R, 

A(x)  is  the  "optimal  upward  directed  rounding",  i.e. 
A(x):=  min  {js.  i £ c M,  £ 2 x)  for  all  x e R • 


(2.5) 


If  p is  an  optimal  rounding,  we  define 

d(T ) :=  ^(t^i  t^,  tR)  :=  (Vt^j  pt^,  Atp)  e (2*6) 
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and  obtain 


Theorem;  If  T is  a triplex  number,  then  4(T)  is  also  a triplex  number. 

Proof;  We  must  show  that  Vti  ptu  <_  AtR.  With  :=  max  { ! x eH,  x £ t^) 
£ tM  we  find  by  (2*4),  (2.3b)  and  (2*3a) 

VtL  = mx  = pm^  <.  ptM  . 

Quite  similarly,  we  have  with  m2  : = min{x  I i tM,  j>  tR)  2.  tM  the  relation 
tj^  Sl  ®2  = ®2  = At  R.  [ ] 

Remark:  A similar  argument  will  show  that  a triplex  number  will  round  to 
a machine  triplex  number  whenever  the  set  of  machine  numbers  into  which  the 
main  value  is  rounded  contains  the  set  of  machine  numbers  into  which  the 
endpoints  are  rounded*  Thus  we  may  use  higher  precision  for  the  main  value 
than  for  the  endpoints  if  we  so  choose;  however,  we  must  not  use  lower  preci- 
sion (it  is  easy  to  construct  a counterexample). 

Thus,  * if  we  wish  to  use  the  centered  form  (as  indicated  above)  to  train 
precision,  we  must  regard  the  main  value  as  having  at  least  the  same  precision 
as  the  error  term,  with  as  many  trailing  zeros  as  necessary,  and  any  nonzero 
digits  beyond  the  nominal  precision  of  the  main  value  must  be  absorbed  into 
the  error  term  prior  to  rounding  the  triplex  number  to  a machine  triplex  num- 
ber. 


In  practice,  we  will  consider  only  triplex  numbers  of  the  form  (2.1),  and 
we  will  use  the  rounding  defined  by  (2.3),  with  p taken  to  be  roundine  to  the 
nearest  machine  number. 


3.  Basic  operations  with  LEVELS*  numbers 

We  have  seen  in  Section  2 that  a triplex  number  may  be  interpreted  as  a 
closed  interval  with  a special  element  of  this  interval  beine  labeled  as  the 
main  value.  Therefore,  the  operations  which  we  now  define  for  triplex  numbers 
are  all  obtained  from  the  correspond ins:  operations  for  intervals  and  an  addi- 
tional operation  for  the  main  value. 


If  • is  one  of  the  four  standard  arithmetic  operators  +,  -,  x,  /, 
T :=(tLi  tM,  tR)  and  W :r(wL,  wM,  wR)  are  triplex  numbers,  then  T • W is 
fined  to  be  the  triple 

T • W :=  (inf  (t  ‘ w | t e IT,  w t IW|,  tM  1 wM, 

sup  (t  * w ! t e IT,  w e IW} ) (provided  0 4 W for  division) 

This  definition  may  be  re-interpreted  as 
I(T  » W)  : (t  • w ! t £ IT,  w e IW)  = IT  » IW,  and  we  have 

* WM  € I(T  * w>»  since  tM  e IT,  wM  e IW, 

with  the  known  interval  operation  IT  • IW  = I(T  • W). 


and 

de- 


(3.1) 
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More  generally,  if  f is  a real-valued  function  of  k real  variables,  which 


is  defined 

and 

continuous  on  IT ^ x IT2  x •»» 

x ITk,  we 

define 

F(T^ , 

T2> 

•••»  :s(inf{f(t1f  t2, 

4k  ' ti  e 

ITil, 

f(t1M»  42M» 

•••»  tkM)« 

(3.2) 

| sup{  f(  t •( , t2,  ..*,  tk|ti  e IT.^})  ( | 

Again,  (3«2)  may  be  interpreted  with  the  known  united  interval  extension  F of 
f as 

KFd,,  T2,  ....  Tk)  = F(IT1f  IT2 ITk)  : = 

{ f ( t i , t2,  ...,  tk)  | ti  e IT.^},  and  we  have 

f(t-jn|,  t2m,  ...,  tk^)  F(IT^,  IT2,  »».,  ITk)  since  t^  e IT^. 

Thus,  (3*1)  and  (3.2)  define  operations  on  triplex  numbers. 

In  both  definitions  there  is,  as  in  the  corresponding  interval  opera- 
tions, the  unpleasant  task  of  computing  the  inf  and  the  sup.  Compactness  ar- 
guments show  that  these  values  may  be  obtained  as  function  values  of  certain 
points  in  IT  x IW  or  IT ^ x IT2  x ...  x ITk,  respectively.  A complete  discus- 
sion for  the  four  arithmetic  operations  and  for  the  elementary  transcendental 
functions  is  given  in  [I7];  several  relational  operators  are  also  discussed 
there. 

1 We  also  need  operators  which  form  triplex  numbers  from  reals  or  vice-ver- 

sa, give  information  about  the  structure  of  the  interval  defined  by  the 
triplex  number,  and  evaluate  certain  logical  functions. 

F<?riing  JOSS  triplex  numbers: 

Compose  is  a function  of  three  real  variables  giving,  in  this  order,  the 
left  endpoint,  the  main  value,  and  the  right  endpoint  of  the  triplex  number. 

Union  is  a function  of  two  triplex  variables  yielding  the  triplex  number 
T defined  by  ti  :=  min  {t«  , t™},  t^  :=  (t^u  + t2M)  / 2,  tp  :=  max  {t1R, 

t2R}»  Thus,  IT  is  the  smallest  interval  containing  IT1  u IT.,. 

Intersect  is  a function  of  two  triplex  numbers  Tj,  T2  yielding  the 

triplex  T with  t^  • " max  {t^,  * " ^tn  {t^R,  t2p},  ^M  * ~ (1^  ^ tp)  ^ 

2,  if  tL  L tp.  For  tL  > tp,  Intersect  is  not  defined.  Thus,  when  it  is  de- 
fined, IT  is  the  largest  interval  in  IT1  n IT2. 

Measures  of  triplex  numbers:  These  are  obtained  as  the  measures  of  the 

corresponding  intervals  as  defined  by  Ris  [14]  and  used  in  [I7].  We  have,  in 

particular,  Absolute  Value  of  T,  Length  of  T,  Half-length  of  T,  Midpoint  of  T, 

Size  of  T,  Suoremum  of  T,  Infimum  of  T,  Magnitude  of  T,  Mignitude  of  T,  Pivot 
of  T,  Sign  of  T,  and  the  distance  between  two  triplex  numbers,  dist(T,  W)  := 

max  { i tL  - wLl,  |tM  - wMl,  itp  - wR  |}. 
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Logical  operators;  These  provide  the  means  to  test  various  properties  of 
triplex  numbers. 

Bad  is  a logical  function  of  one  triplex  number.  Its  value  is  true  if 
and  only  if  T is  an  improper  triplex  number,  that  is,  if  tL  > tM  or  tp  < tM» 

OK  is  the  logical  inverse  of  Bad;  i.e.,  its  value  is  true  if  and  only  if 
T is  a proper  triplex  number  — that  is,  if  tL  <.  tM  <.  tR. 

Element  o£  is  a logical  function  of  a real  number  x and  a triplex  number 
T.  Its  value  is  true  if  and  only  if  x is  an  element  of  IT. 

Subset  (Superset ) is  a logical  function  of  two  triplex  numbers  T^  and  T2» 
Its  value  is  true  if  and  only  if  IT^  is  contained  in  IT2  (IT2  is  contained  in 
ITt). 


One  main  point  in  implementing  the  definitions  stated  here  is  to  take  ac- 
count of  the  effects  of  rounding.  In  Section  2,  we  have  already  pointed  out 
that  this  rounding  has  to  be  done  so  as  to  ensure  that  the  left  endpoints  are 
rounded  downward,  the  right  endpoints  are  rounded  upward,  and  the  main  values 
are  rounded  to  the  nearest  machine  number.  If  F is  a triplex  function  as  de- 
fined in  (3*2),  therefore,  we  obtain  the  rounded  value  as 


E(T-p  T2,  ...,  Tk)  ;=  t5(F(T^,  T2,  ...,  Tk)  ; = 
(7inf{f(t^«  t2,  ...,  tk)  i t^  c ITjJ, 
pf(t1M,  t2M,  ...,  tkM), 

As  up  ( f ( 1 1 » 1 2 , • f tk)  ! ti  6 IT.}) 


(3.3) 


where  V denotes  the  optimal  downward  directed  rounding,  A denotes  the  op- 
timal upward  directed  rounding,  and  p denotes  the  optimal  rounding  to  the 
nearest  machine  number. 


In  '3.3)  it  is  possible  to  evaluate  the  endpoints  as  for  ordinary  inter- 
vals (see  for  example  [1?])  and  the  main  values  may  be  computed  directly,  pro- 
vided that  (a)  the  computation  is  carried  out  with  at  least  as  much  precision 
as  is  used  for  the  endpoints,  and  (b)  the  real  arithmetic  used  for  the  compu- 
tation rounds  the  result  to  the  machine  number  nearest  the  true  result  (many 
computers  do  not  do  this). 

Another  problem  is  that  of  out-of-range  numbers.  We  refer  the  reader  to 
the  discussion  given  in  Cl7].  The  only  additional  difficulty  arising  in 
triplex  numbers  is  the  inclusion  of  the  main  value  tM»  Since  everv  improperly 
computed  triplex  number  should  lead  to  an  error  message,  we  must  check  each  of 
tLi  t^,  and  tR  to  see  whether  it  was  computed  successfully  or  whether  one  of 
the  errors  overflow,  infinity,  or  underflow  was  encountered.  The  possible  er- 
rors are  documented  in  Appendix  3. 


4.  Using  the  TRIPLEX  arithmetic  package: 

This  section  is  designed  to  be  the  user’s  guide  for  the  triplex  arithme- 
tic package.  The  package  is  most  easily  used  with  the  aid  of  the  AUGMENT  pre- 
compiler;  this  use  will  be  described  first*  Following  that,  we  discuss  the 
use  of  the  package  without  benefit  of  the  precompiler;  essentiall-:,  this  is 
just  a discussion  of  how  the  user  should  go  about  doing  the  things  that  the 
precompiler  would  do  were  it  available* 

4. 1 Using  the  package  with  the  AUGMENT  precompiler: 

Declaring  TRIPLEX  variables:  If  X is  a TRIPLEX  variable,  Y is  a TRIPLEX 
vector  of  length  10,  and  Z is  a 10  x 10  TRIPLEX  matrix,  the  following  type 
declaration  will  both  establish  them  as  TRIPLEX  variables  and  reserve  storage 
for  them: 

TRIPLEX  X,  Y(10) , Z( 1 0, 10) 

Restrictions: 

Identifiers  should  not  begin  with  the  prefixes  "TPX'  or  "BPA';  this 
will  help  avoid  conflicts  with  the  names  of  subroutines  in  the  package. 

If  the  FORTRAN  compiler  limits  the  number  of  dimensions  of  an  arrav, 
that  limit  must  be  decreased  by  1 for  TRIPLEX  variables;  AUGMENT  will  de- 
clare an  undimensioned  TRIPLEX  variable  as  a vector. 

Options:  None 

Assigning  constant  values  io  TRIPLEX  variables : This  may  be  accomplished 
by  a statement  which  assigns  the  value  of  an  appropriate  Hollerith  string  to 
the  TRIPLEX  variable.  For  example,  to  assign  the  value  (.1,  .1,  .1)  to  the 
TRIPLEX  variable  X,  one  would  write 

X = 1 3H ( . 1 , .1,  .1)$ 


Restrictions: 

The  Hollerith  string  must  represent  a TRIPLEX  as  defined  in  Appendix 
1,  except  that  embedded  blanks  are  permitted  and  will  be  ignored. 

The  Hollerith  string  must  be  terminated  by  a field  separator:  $,  #, 
or  = . 

The  length  of  the  Hollerith  string  is  limited  to  the  length  of  the 
formatted  read  buffer  (132  characters  in  the  UNIVAC  1110  version). 

The  number  of  significant  digits  permitted  for  each  endpoint  of  the 
triplex  number  is  limited  to  the  value  of  ESDMAX  (60  in  the  UNIVAC  1110 
implementat ion) • 


The  user  should  be  aware  that  a statement  of  the  form 


( 

L 


X 


1 


will  cause  a value  to  be  assigned  to  X , but  the  resulting  interval  of 
the  triplex  number  may  not  contain  .1.  The  above  statement  will  cause  X 
to  be  set  equal  to  a degenerate  triplex  number  each  of  whose  endpoints  is 
the  REAL  approximation  to  .1;  in  order  to  obtain  a rigorous  triplex  num- 
ber approximation  to  .1  on  most  computers,  the  triplex  number  would  need 
to  be  nondegenerate*  Thus  shortcutting  the  conversion  from  Hollerith  to 
triplex  can  be  dangerous. 

Options: 

Quoted  Hollerith  literals  may  be  used  if  the  host  compiler  accepts 
them* 


If  the  host  compiler  generates  a sentinel  for  Hollerith  literals,  and 
if  the  TPXUPK  primitive  recognizes  this  sentinel  (as  is  the  case  in  the 
UNIVAC- 1110  implementation),  the  terminal  field  separator  may  be  omitted. 

Implicit  conversations  from  Hollerith  to  TRIPLEX  are  allowed  (e*g.,  a 
Hollerith  literal  may  appear  in  an  arithmetic  expression). 

Reading  TRIPLEX  variables  --  free  format : The  free  format  read  will  ob- 
tain the  next  data  field  from  the  input  stream  on  the  specified  unit,  convert 
it,  and  store  the  resulting  value  in  the  specified  triplex  variable.  This  may 
be  accomplished  by  a statement  of  the  form 

CALL  TPXRDF(UNIT,X) 

Restrictions: 

Only  one  value  is  read  by  each  call  to  TPXRDF. 

The  basic  package  recognizes  units  5 (standard  input)  and  0 (standard 
reread ) * 

The  Hollerith  string  read  by  TPXRDF  must  represent  a triplex  number 
as  defined  in  Appendix  1. 

The  length  of  the  Hollerith  string  is  limited  to  the  length  of  the 
formatted  read  buffer  (132  characters  in  the  UNIVAC  1110  version). 

The  number  of  significant  digits  permitted  for  each  endpoint  of  the 
triplex  number  is  limited  to  the  value  of  ESDMAX  (60  in  the  UNIVAC  1110 
version) . 

The  end  of  an  input  record  does  not  terminate  the  input  stream. 

A call  to  the  free-format  read  routine  with  a different  unit  number 
specified,  or  any  call  to  the  formatted  read  routine  (see  below),  will 
terminate  the  input  stream  on  the  current  unit.  Once  the  incut  stream 
has  been  terminated,  TPXRDF  begins  a new  input  stream  with  a new  record. 
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TPXRDF  will  read  only  from  units  designated  as  read  units  bv  the  in- 
formation in  UNITBL  (see  "Options"  below). 


Options: 

The  standard  unit  designations  may  be  changed  by  altering  UNITBL, 
which  is  located  in  the  COMMON  block  TPXCCM.  (This  change  would  normally 
be  accomplished  at  the  time  of  adaptation,  however.) 

Descriptions  of  additional  I/O  may  be  entered  in  UNITBL.  To  do  this 
NUNITS  must  be  changed  to  allow  the  I/O  routines  to  scan  all  information 
in  the  table  (increase  NUNITS  by  1 for  each  unit  added).  Each  unit  de- 
scription consists  of  a row  of  UNITBL.  Column  1 is  the  unit  number;  Col- 
umn 2 is  the  length  of  the  record  (limited  to  the  length  of  the  free  for- 
mat read  buffer,  which  is  132  characters  in  the  1110  version);  Column  3 
is  the  number  of  characters  in  each  record  to  be  scanned;  Column  4 is  0 
if  the  record  does  not  contain  a carriage  control  character  and  1 if  it 
does;  Column  5 is  1 if  the  unit  is  read  only,  2 if  write  only,  and  3 if 
read/write. 

Units  may  be  referenced  by  their  Dosition  in  the  table  rather  than 
their  logical  unit  numbers:  -1  refers  to  the  unit  described  in  the  first 
row  of  the  table,  and  so  on. 

The  number  of  characters  of  each  record  which  are  actually  scanned 
may  be  changed  if  desired;  this  allows  one  to  use  only  the  first  72  char- 
acters in  a card  image,  for  example.  Change  UNITBL(1,3)  to  the  aporoori- 
ate  value  to  accomplish  this. 

The  end  of  a record  can  be  treated  as  the  end  of  the  data  field  by 
setting  UNITBL(1,  3)  to  one  greater  than  the  length  of  the  physical  rec- 
ord and  inserting  a field  termination  character  in  the  corresDonding 
character  of  IBUFRE.  The  length  of  this  buffer  must  not  be  exceeded  in 
doing  this,  however. 

The  remainder  of  the  current  input  record  may  be  skiDDed  bv  setting 
OLDUNT  to  -1  prior  to  calling  TPXRDF. 

Each  input  record  may  be  echoed  on  the  standard  print  unit  as  it  is 
read.  To  accomplish  this,  set  ECHOS  to  .TRUE,  in  COMMON  block  TPXCRD. 

Each  data  field  may  be  echoed  on  the  standard  print  unit  as  it  is 
converted.  To  do  this,  set  ECHOD  to  .TRUE,  in  common  block  TPXCRD. 

Delimiters  for  triplex  data  may  be  altered  as  desired  by  changing 
the  appropriate  value  in  ICHR,  located  in  COMMON  block  TPXCCM.  Note  that 
this  will  also  affect  the  Hollerith  strings  used  to  assign  constant  val- 
ues to  triplex  variables.  If  a character  is  duplicated  in  this  arrav, 
its  interpretation  will  be  governed  by  its  first  appearance. 


Data  fields  may  be  separated  by  blanks  (as  many  as  desired).  Blanks 
occurring  between  matching  parenteses  will  be  ignored. 


Commas  within  matching  pairs  of  parentheses  are  regarded  as  endpoint 
separators  (unless  the  ICHR  array  is  redefined). 

A field  may  be  terminated  by 

1*  A blank  which  is  neither  a leading  blank  nor  located  be- 
tween matching  parentheses; 

2.  Any  of  the  characters  '#',  or 

3.  A comma  occurring  outside  of  matching  Darentheses; 

4.  Any  nonblank  character  following  a matching  right  parenthe- 
sis; if  this  is  not  one  of  the  characters  mentioned  in  (2) 
or  (3)  above,  then  it  will  be  regarded  as  the  first  charac- 
ter of  the  next  field. 

Any  null  field  or  subfield  is  taken  to  represent  the  number  0. 

A field  containing  a single  number  will  be  converted  to  the  smallest 
triplex  number  containing  that  number.  The  resulting  triplex  number  may 
or  may  not  be  degenerate. 

If  a field  contains  only  two  numbers,  these  are  ta‘  ^.n  to  be  the 
endpoints  of  the  triplex  numbers  and  the  main  value  is  taken  to  be  the  mid- 
point of  the  interval  thus  specified. 

Reading  TRIPLEX  variables  — formatted:  The  formatted  read  routine  reads 
and  converts  a vector  of  data;  the  vector  may,  of  course,  be  of  length  1*  The 
user  supplies  a format,  which  is  an  array  of  length  3:  The  first  element  is 
the  number  of  data  fields  in  each  record;  the  second  element  is  the  number  of 
characters  to  be  skipped  before  beginning  the  scan  on  each  field;  and  the 
third  is  the  number  of  characters  in  each  data  field.  The  calling  seouence  is 

CALL  TPXRD  (UNIT,  FMT,  A,  N) 

where  UNIT  is  the  unit  number,  FMT  is  the  format  as  described  above;  A is  the 
first  location  of  the  vector  into  which  the  data  is  to  be  read;  and  N is  the 
length  of  the  vector. 

Restrictions: 

The  basic  package  recognizes  units  5 (standard  input)  and  0 (re- 
read ) • 

The  contents  of  each  field  must  represent  a triplex  number  as  de- 
fined in  Appendix  1.  Field  termination  characters  are  not  permitted. 

The  length  of  each  record  is  limited  to  the  length  of  the  formatted 

read  buffer. 

The  number  of  significant  digits  permitted  for  each  endpoint  is  lim- 
ited to  the  value  of  ESDMAX. 
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TPXRD  will  read  only  from  units  designated  as  read  units  in  UNITBL. 
Options: 

The  standard  unit  designations  may  be  changed  as  in  TPXRDF. 

I/O  units  may  be  added  to  UNITBL;  see  discussion  of  TPXRDF. 

Units  may  be  referenced  by  their  positions  in  UNITBL  as  discussed 
under  TPXRDF. 

The  number  of  characters  in  a record  which  are  actually  used  may  be 
altered  as  discussed  under  TPXRDF* 

Input  data  may  be  echoed  as  discussed  under  TPXRDF* 

Delimiters  for  triplex  data  may  be  altered  as  discussed  under 
TPXRDF* 

Note:  Any  of  the  above  modifications  will  affect  both  of  the 
read  routines*  Blanks  may  be  embedded  in  the  field;  they  will 
be  ignored. 

The  use  of  parentheses  is  optional* 

A comma  will  be  interpreted  as  an  endpoint  separator  regardless  of 
whether  parentheses  are  used* 

If  information  is  read  from  a unit  on  which  a carriage  control  char- 
acter is  indicated,  the  first  character  of  the  record  will  be  ignored. 

Computing  with  TRIPLEX  variables:  Expressions  involving  TRIPLEX  varia- 
bles are  written  in  standard  FORTRAN  syntax,  just  as  though  TRIPLEX  were  a 
standard  FORTRAN  data  type*  A list  of  the  operations  and  functions  available 
in  this  package  may  be  found  in  Appendix  2. 

The  list  of  operations  and  functions  available  includes  all  appropriate 
ANSI  standard  FORTRAN  operations  and  functions  — i.  e.,  all  which  have  natu- 
ral triplex  extensions.  Other  operators  and  functions  peculiar  to  triplex 
arithmetic  and  extensions  to  triplex  numbers  are  implemented;  examples  include 
the  union  of  two  intervals,  the  midpoint  of  an  interval,  and  various  measures 
of  the  size  of  an  interval.  Relational  operators  are  also  implemented,  but 
these  take  on  different  meanings  in  the  context  of  triplex  arithmetic,  and 
therefore  have  been  given  slightly  different  operator  symbols* 

Restrictions: 

Variable  and  subprogram  names  beginning  with  TPX  or  BPA  should  be 
avoided  to  preclude  conflicts  with  the  package. 

If  COMMON  block  declarations  are  included  in  the  program  for  the 
purpose  of  exercising  various  options,  then  the  variable  names  associated 
with  these  COMMON  blocks  must  likewise  be  avoided. 
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Options: 

Mixed  mode  expressions  are  permitted,  although  their  use  is  discour- 
aged due  to  the  high  probability  of  introducing  hidden  errors*  For  exam- 
ple, the  expression 

Y = 0.1  • X 

where  X and  Y are  TRIPLEX  variables,  will  not  yield  a correct  value 
for  Y:  0*1  will  first  be  converted  to  REAL  by  the  compiler,  and  AUGMENT 
will  then  cause  that  REAL  number  to  be  converted  to  a degenerate  triplex 
number  which  will  not  contain  *1  (unless  the  computer  base  is 
decimal -related) ; multiplication  will  then  take  place  using  the  erroneous 
triplex  number*  If  mixed  mode  expressions  are  used,  TRIPLEX  takes  pre- 
cedence over  REAL,  INTEGER,  and  DOUBLE  PRECISION  (and  MULTIPLE,  for  the 
multiple  precision  version  of  the  package)  data  types*  The  data  type  to 
be  used  in  evaluating  each  subexpression  is  determined  by  normal 
heierarchy  rules. 

New  special  function  routines  may  be  added.  If  this  is  done,  the 
following  conventions  and  rules  should  be  observed  in  order  to  preserve 
the  integrity  of  the  package: 

1*  The  name  of  the  new  routine  should  begin  with  TPX.  Up  to  three 
more  characters  may  be  appended  to  form  the  name  of  the  routine. 
Conflicts  with  existing  names  must  be  avoided. 

2.  The  calling  sequence  for  the  new  routine  should  adhere  to  the 
standards  of  the  package:  the  arguments  are  listed  first,  fol- 
lowed by  the  result.  Every  effort  should  be  made  to  minimize 
the  amount  of  information  passed. 

3.  It  is  the  user's  responsibility  to  insure  that  the  endpoints  of 
the  resulting  triplex  numbers  are  computed  and  rounded  correct- 
ly. All  existing  package  routines  are  available  to  the  user  in 
writing  the  new  routine,  just  as  they  are  in  writing  an  applica- 
tions program. 

4.  If  the  new  routine  is  to  be  invoked  by  the  AUGMENT  precompiler, 
the  description  deck  for  AUGMENT  must  be  modified  to  include  the 
new  routine.  Consult  [5]  for  details. 

5.  If  errors  are  possible  in  the  new  routine,  TPXRAP  should  be  used 
to  handle  the  error  message  rather  than  including  an  error 
print-  out  in  the  new  subprogram.  This  is  done  as  follows: 

a.  The  suffix  of  the  name  of  the  new  routine  is  stored  in  the 
next  available  location  in  the  NAMES  array. 

b.  Appropriate  values  to  indicate  the  data  types  of  the  argu- 
ments to  the  new  routine  are  stored  in  TYPA,  TYPB,  and 
TYPR.  Note  that  it  is  assumed  that  the  routine  will  have 
at  most  two  arguments  and  one  result. 
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c.  The  TPXCOM  declarations  should  be  included  in  the  new  rou- 
tine. If  any  argument  or  result  is  of  a tyDe  other  than  a 
standard  type,  BPA,  TRIPLEX  or  EXTENDED,  a variable  of  the 
appropriate  type  must  be  defined  and  EQUIVALENCEd  to  TA, 
TB,  or  TR  respectively;  storage  of  the  data  in  the  COMMON 
block  is  then  done  by  using  the  variable  name  of  the  aporo- 
priate  type. 

d.  TPXRAP  should  be  called  with  the  statement 

CALL  TPXRAP 

immediately  prior  to  the  RETURN  statement. 

e.  Immediately  prior  to  the  call  on  TPXRAP,  the  values  of  the 
arguments  must  be  stored  in  TA  and  TB  and  the  value  of  the 
result  must  be  stored  in  TR.  TPXFLT  must  be  set  to  the  ap- 
propriate fault  indicator  (If  no  existinng  fault  value  is 
approrpirate,  a new  fault  value  will  need  to  be  defined  and 
TPXRAP  will  have  to  be  modified  to  provide  the  proper 
printout  — this  constitutes  a major  modification  of  the 
package).  Finally,  the  value  of  ID  must  be  set  equal  to 
the  index  of  the  name  of  the  new  routine  in  the  NAMES  ar- 
ray. 

Writing  TRIPLEX  variables:  The  write  routine  will  convert  a vector  (pos- 
sibly of  length  1)  of  triplex  variables  to  external  format  and  write  it  on  the 
specified  unit  according  to  the  given  format.  The  format  is  now  an  integer 
array  of  length  4:  the  first  three  elements  of  the  array  are  the  same  as  for 
the  formatted  read,  except  that  unused  characters  between  fields  are  filled 
with  blanks;  FMT(4)  is  a Hollerith  carriage  control  character  ('  'or  'O'  in 
the  1110  implementation)  for  use  where  approprite.  The  calling  sequence  is 

CALL  TPXWR  (UNIT,  FMT , A,  N) 

where  the  arguments  are  as  described  in  the  formatted  read. 

The  external  representation  of  each  triplex  number  is  guaranteed  to  con- 
tain the  triplex  number,  and  is  the  triplex  number  with  the  smallest  interval 
representable  in  the  given  format  which  does  so. 

Restrictions: 

The  basic  package  recognizes  units  6 (standard  print  unit)  and  1 
(standard  punch  unit). 

The  length  of  each  record  is  limited  to  the  length  of  the  write  buf- 
fer (132  characters  in  the  UNIVAC  1110  implementation). 

The  width  of  each  data  field  specified  by  the  format  must  be  at 
least  great  enough  to  permit  the  package  to  convert  one  significant  dig- 
it. In  the  1110  single  precision  version,  this  is  15  digits,  assuming  a 
2-digit  exponent.  Add  two  characters  for  each  digit  of  exponent  above  2. 


- 14  - 


I 


The  number  of  significant  digits  allowed  for  each  endpoint  is  limit- 
ed to  the  value  of  ESDMAX  (60  in  the  1110  implementation)* 

TPXWR  will  write  only  on  units  designated  as  write  units  in  UNITBL* 

Note:  if  any  of  the  restrictions  are  violated,  TPXWR  will  use 
a standard  format  and/or  the  standard  print  unit,  as  necessary, 
to  insure  that  the  designated  information  is  written  out  in 
some  form* 


Options: 

The  standard  unit  designations  may  be  changed  (see  TPXRDF)* 

I/O  units  may  be  added  to  UNITBL  (see  TPXRDF)* 

Units  may  be  referenced  by  their  positions  in  UNITBL  (see  TPXRDF)* 

Delimiters  for  triplex  data  may  be  changed  by  altering  the  appropri- 
ate values  in  OCHR,  located  in  COMMON  block  TPXCCM* 

Carriage  control  specifications  may  be  changed  by  altering  the  ap- 
propriate element  of  UNITBL* 

Errors : The  package  is  designed  to  detect  all  errors  as  they  occur*  The 
TPXRAP  routine  will  print  an  error  message  and  halt  the  computation  except  in 
those  cases  where  a viable  alternative  exists*  Those  cases  are  few  indeed; 
they  comprise  arithmetic  underflows  (where  the  offending  value  is  set  eaual  to 
zero  or  to  the  properly  signed  number  of  smallest  magnitude,  as  aoprororiate) 
and  errors  occurring  on  output  (where  the  write  routine  uses  standard  modes  of 
output  in  lieu  of  the  erroneous  information)*  In  the  former  case,  computation 
proceeds  without  notice  to  the  user;  in  the  latter  case,  a message  is  produced 
on  the  standard  print  unit  after  output  is  complete* 

The  possible  errors  and  the  response  to  each  are  listed  in  Appendix  3* 

Restrictions:  None. 

Options: 

The  default  response  to  any  fault  except  80  (Conversion  arrav  over- 
flow) may  be  changed*  To  alter  the  response  to  Fault  number  I,  change 
the  value  of  M0N(I+1)  as  desired*  This  array  is  located  in  TPXRCM* 
MON(80)  should  never  be  changed,  since  this  could  result  in  looping  under 
some  circumstances* 

TPXRAP  may  be  used  as  a trace  routine  of  sorts  by  changing  the  re- 
sponse to  a successful  operation  [M0N(1)3*  This  will  cause  TPXRAP  to 
produce  output  even  after  a successful  operation;  however,  this  output 
is,  of  course,  limited  to  those  routines  which  call  TPXRAP*  In  cases 
where  errors  are  not  possible,  TPXRAP  is  usually  bypassed* 
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Producing  an  object  program;  The  generation  of  an  object  program  is  a 
two-step  process* 

1*  Use  AUQ4ENT  to  translate  the  source  program  into  a FORTRAN  program 
compatible  with  the  host  FORTRAN  compiler*  This  can  be  accomplished 
with  a run  stream  of  the  following  type  (see  Appendix  6): 

invoke  AUGMENT 

[DESCRIPTION  decks  for  TRIPLEX  and  BPA  (see  Appendix  5)] 

•BEGIN 

[source  program] 

•END 

AUGMENT  will  write  the  translated  program  on  unit  20*  (Note:  AUG- 

MENT also  uses  Unit  21,  for  scratch  work*) 

2*  Compile  the  output  of  AUGMENT  using  the  standard  FORTRAN  compiler; 
then  load  and  execute  as  with  any  FORTRAN  program*  The  user  must  in- 
sure that  the  TRIPLEX  package  library  is  available  to  the  linkage  ed- 
itor, and  that  the  BLOCK  DATA  modules  are  included  with  the  program 
by  the  linkage  editor* 

The  run  stream  for  the  UNIVAC  1110  version  at  the  University  of  Wisconsin 
is  shown  in  Appendix  6. 


4.2  uai&g  Basta&a  .vtiUisui  tens  fit  sL  augment: 

If  necessary,  the  TRIPLEX  package  may  be  used  directly,  without  using 
AUGMENT  to  preprocess  the  source  code*  If  this  is  done,  certain  of  the  things 
that  AUGMENT  would  normally  take  care  of  must  be  done  instead  bv  the  user*  In 
this  section,  we  indicate  the  changes  that  are  necessary  in  the  previous  in- 
structions if  AUGMENT  is  not  available*  These  remarks  apply  primarily  to  the 
single  precision  version  of  the  triplex  package*  Suitable  modification  of 
these  instructions  will  render  them  applicable  to  the  multiple  precision  ver- 
sion; however,  due  to  the  relative  complexity  of  the  multiple  precision  ver- 
sion, its  use  in  the  absence  of  AUGMENT  is  discouraged* 

Declaring  TRIPLEX  variables:  The  number  of  words  of  storage  that  must  be 
reserved  for  a triplex  variable  will  be  three  times  the  number  of  words  needed 
to  store  each  endpoint*  In  the  UNIVAC  1110  implementation,  each  triplex  vari- 
able requires  three  words  of  storage*  For  a particular  installation  or  aDDli- 
cation,  however,  the  package  may  have  been  revised  to  produce  greater  accura- 
cy, so  the  local  documentation  should  be  consulted* 

The  package  assimes  that  the  words  containing  a triplex  number  will  occu- 
py consecutive  locations  in  storage* 


I 


The  standard  data  type  used  to  reserve  storage  for  TRIPLEX  variables  is 
irrelevant  so  long  as  consistency  is  maintained.  AUGMENT  has  been  instructed 
to  declare  TRIPLEX  variables  as  REAL  arrays  of  length  3*  If  the  latter  con- 
vention is  used,  a dimensioned  array  of  TRIPLEX  variables  must  be  declared  as 
a REAL  array  of  one  more  dimension,  the  first  dimension  alwavs  being  3* 


If  X is  a TRIPLEX  variable,  Y is  a TRIPLEX  vector  of  length  10,  and 
is  a 10  x 20  TRIPLEX  matrix,  they  could  be  declared  in  the  following  wavs 


REAL  X(3) , Y(3,  10),  Z(3,  10,  20) 


Assigning  constant  yalUfeS  12  TRIjPkfiX  variables:  This  must  now  be  accom- 
plished by  a call  on  the  subroutine  TPXASG.  Depending  on  the  data  type  used 
to  declare  the  triplex  variables  and  the  whims  of  the  compiler,  it  may  be  nec- 
essary to  refer  to,  say,  X(1)  rather  than  X.  The  rules  for  forming  the  Hol- 
lerith strings  are  the  same  whether  or  not  AUGMENT  is  used. 


If  TRIPLEtXes  are  declared  as  REAL  arrays  of  length  3,  then  in  order  to 
assign  the  value  of  (.1,  .1,  .1)  to  X,  one  would  write: 

CALL  TPXASG  ( 1 3H ( . 1 , .1,  .1)$,  X(D) 


Restrictions:  No  change. 


Options:  No  change. 


Reading  TRIPLEX  variables  — free  format : No  change,  unless  the  compiler 
requires  explicit  subscripts  for  dimensioned  arrays  and  TRIPLEX  variables  are 
declared  as  arrays.  In  that  case,  the  subscripts  must  be  furnished. 


Reading  TRIPLEX  variables  — formatted:  See  comments  for  free  format 

read. 

Computing  with  TRIPLEX  variables:  Here,  several  pitfalls  await  the 
unwary  user.  First,  of  course,  is  the  necessity  of  parsing  the  arithmetic  ex- 
pressions. While  this  is  of  vital  importance,  we  assume  that  the  user  knows 
how  to  do  it  correctly. 

Next  is  the  necessity  of  expressing  triplex  variables  in  a form  compati- 
ble with  the  declarations  and  the  restrictions  of  the  compiler.  We  have 
touched  on  this  above,  and  will  not  iterate  here. 

Third  is  the  problem  of  making  sure  that  we  are  doing  exactly  what  we 
want  to  do.  For  example,  if  we  wish  to  multiply  the  TRIPLEX  variable  X bv 
2.0,  we  can  not  write 

CALL  TPXMUL  (X,  2.0,  Y) 

since  TPXMUL  could  multiply  X by  a "triplex"  whose  left  endpoint  was  2.0  and 
whose  main  value  and  right  endpoint  were  whatever  haopended  to  bv  lying  around 
in  the  next  two  storage  locations.  The  package,  of  course,  has  no  way  of  de- 
tecting such  blunders;  the  user  alone  is  responsible  for  avoiding  them. 
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Thus  the  user  is  responsible  for  employing  explicit  data  type  conversions 
where  they  are  needed,  and  insuring  that  all  evaluations  of  expressions  are 
performed  using  the  correct  data  types. 

Restrictions:  No  change 

Options:  No  change. 

MrlUnfi  .TRIPLEX  v^ria^:  See  comments  for  free  format  read. 

Errors:  No  change. 

Producing  an  object  program:  The  AUGMENT  phase  will,  of  course,  not  ap- 
ply. The  remaing  remarks  are  valid. 


4. 3 Using  BP A variables  in  the  source  program: 

The  data  type  of  the  endpoints  of  a triplex,  and  of  the  results  of  cer- 
tain triplex  functions,  is  a nonstandard  type  named  BPA.  If  BPA  format  is  the 
same  as  REAL,  computations  with  these  quantities  can  usually  be  executed  using 
REAL  arithmetic.  If  a REAL  variable  is  set  equal  to  a BPA  expression,  AUGMENT 
will  perform  an  automatic  type  conversion;  if  AUGMENT  is  not  used,  BPACBR 
should  be  called  to  effect  this  conversion.  If  mixed  mode  expressions  involv- 
ing BPA  and  standard  data  types  are  written,  AUGMENT  will  cause  the  expression 
to  be  evalued  using  BPA  arithmetic,  converting  the  result  to  a standard  data 
type  if  necessary. 

Computation  using  BPA  routines  may  be  included  in  the  source  program,  ei- 
ther implicitly,  as  above,  or  explicitly.  THE  USER  IS  CAUTIONED  THAT  MOST  BPA 
ROUTINES  REQUIRE  THAT  A ROUNDING  OPTION  BE  STORED  IN  BPACOM  PRIOR  TO  INVOKING 
THE  ROUTINES  (see  [I7]).  Thus,  for  example,  in  order  to  add  BPA  numbers  A and 
B,  using  upward-directed  rounding,  and  store  the  result  in  C,  one  would  need 
the  following  statements  in  the  program: 

COMMON  /BPACOM/  ...  (declarations  from  Appendix  4) 

BPA  A,  B,  C 


OPTION  = RDU 
C s A ♦ B 

If  OPTION  is  not  set  explicitly,  it  will  have  whatever  value  it  may  pre- 
viously have  had,  and  the  rounding  may  not  be  as  expected.  This  could  result 
in  erroneous  results. 


5.  Design . implementation,  ao£  adaptation  other  hardware: 

The  current  triplex  arithmetic  package  was  a straight  forward  modification 
of  the  interval  arithmetic  package  discussed  in  [ 1 7 ] • In  light  of  the 
foregoing  discussion,  it  should  not  be  surprising  that  all  that  was  necessary 
was  to  add  code  to  all  of  the  interval  routines  to  compute  the  main  values, 
adjust  formats  to  accommodate  the  third  component,  change  the  representation 
of  intervals  to  that  of  triplexes,  and  add  a few  routines  to  handle  the  main 
values* 

Since  the  interval  arithmetic  package  had  been  written  with  the  aid  of 
AUGMENT,  using  the  nonstandard  data  type  INTERVAL,  modification  of  the  package 
was  not  at  all  difficult*  We  first  added  a field  function  MAIN,  which  either 
inserts  or  extracts  the  main  value  of  the  triplex  number  according  as  it  ap- 
pears on  the  left  or  right  side  of  the  = sign*  This  is  known  as  a field  func- 
tion in  the  jargon  of  AUGMENT,  and  is  permitted  by  AUGMENT  whether  or  not  such 
things  are  allowed  by  the  FORTRAN  compiler*  The  function  MAIN  was  entirely 
analogous  to  the  field  functions  INF  and  SUP,  which  we  had  used  in  the 
inplementation  of  the  interval  arithmetic  package  (see  Appendix  2). 

Having  added  the  MAIN  function,  it  was  mechanical  to  add  the  lines  of 
code  to  each  routine  to  compute  the  main  value  of  the  triplex  number*  For  ex- 
ample, in  the  square  root  routine,  the  computation  of  the  main  value  looks 
like 

MAIN(R)  = SQRT(MAIN(A) ) 

As  noted  above,  errors  may  occur  in  computation  of  the  main  values  as  well  as 
in  the  computation  of  the  endpoints  of  triplex  numbers.  We  therefore  revised 
the  table  of  faults  to  accommodate  the  additional  possibilities*  The  revised 
table  is  shown  in  Appendix  3*  This  also  necessitated  some  chances  to  the  er- 
ror-handling routine  INTRAP,  of  course* 

Next,  we  renamed  the  modules  of  the  interval  package  by  suitable  applica- 
tion of  a text  editing  program;  we  merely  replaced  the  prefix  INT  with  the 
prefix  TPX  on  each  module*  Again,  since  we  were  usinc  AUGMENT  to  venerate 
calls  on  other  package  modules,  there  was  a minimum  of  work  to  do  in  makine 
these  changes  — the  most  pervasive  change  was  the  call  to  the  error-  handling 
routine,  which  appeared  in  nearly  every  routine* 

Changing  the  representation  of  intervals  to  that  of  triplex  numbers  was 
very  easy*  In  each  routine,  interval  data  had  been  declared  type  INTERVAL. 
We  simply  changed  INTERVAL  to  TRIPLEX  in  the  package  modules,  then  supplied 
AUGMENT  with  a new  description  deck  which  defined  the  structure  of  triplex 
numbers. 

We  did,  of  course,  need  to  code  the  functions  which  implemented  the 
field  function  MAIN,  and  revise  the  primitive  functions  which  had  defined  in- 
tervals in  terms  of  arrays  of  BPA  numbers  (such  as  INF  and  SUP)*  These  were 
not  difficult  tasks* 
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Perhaps  the  most  difficult  of  all  of  the  tasks  were  revision  of  the  input 
routine  to  accommodate  triplex  numbers  and  alteration  of  the  output  formats. 
Even  so,  these  modifications  took  less  than  one  day. 

Obviously,  then,  the  triplex  package  described  here  is  very  closely  re- 
lated to  the  interval  arithmetic  package  described  in  [I7].  The  design  phi- 
losophy of  the  triplex  package  is  identical  to  that  of  the  interval  package, 
and  need  not  be  discussed  here.  The  underlying  package  for  BPA  arithmetic  is 
the  same  package  in  all  particulars,  so  we  refer  the  reader  interested  in  that 
part  of  the  triplex  package  directly  to  [I7]* 

Nearly  all  of  Section  5 of  [I7]  (Design  and  implementation  of  the  [inter- 
val] package)  and  Section  6 (adapting  the  [interval]  package  to  other  hard- 
ware) transfers  directly  to  the  triplex  package  if  the  reader  mentally  makes 
the  same  changes  in  the  text  that  we  made  in  the  program:  interval  to 
triplex;  the  prefix  INT  to  the  prefix  TPX  in  names  of  package  routines  (e.g., 
INTCOS  would  become  TPXCOS);  definition  of  triplex  numbers  as  BPA  arrays  of 
length  3 instead  of  interval  numbers  as  BPA  arrays  of  length  2.  Of  course, 
there  are  other,  less  mechanical  changes  (such  as  the  definition  of  -T,  where 
T is  a triplex  number:  if  T = (a,  b,  c),  then  -T  :=  (-c,  -b,  -a)  )•  The  flow 
chart  of  a typical  triplex  function  would  need  to  include  steps  for  calcula- 
tion of  the  main  value.  And  the  state  table  for  the  free-format  read  would 
need  to  be  revised  to  admit  the  possibility  of  triples  instead  of  pairs.  But 
none  of  these  things  are  very  subtle,  and  the  reader  should  not  have  to  suffer 
the  indignity  of  having  them  spelled  out  in  detail. 


6.  The  multiple 


oq.  version  £f  the  J.rip.le^  Paskage: 


In  the  foregoing  discussion,  we  have  considered  primarily  the  single  pre- 
cision version  of  the  Triplex  Package  (although  the  instructions  given  for  use 
of  the  package  with  AUGMENT  are  equally  applicable  to  the  multiple  precision 
version).  The  single  precision  version  was  initially  developed  for  the  UNIVAC 
1110  computer,  but  can  be  transferred  to  other  hardware  without  a great  deal 
of  reprogramming.  Indeed,  most  of  the  necessary  work  is  in  the  BPA  portion  of 
the  package,  and  this  has  already  been  adapted  to  several  other  machines.  See 
[ 18]  for  details. 

The  multiple  precision  arithmetic  package  of  Brent  [3]  is  a completely 
portable,  arbitrary-precision  package  which  includes  modules  to  evaluate  all 
of  the  standard  mathematical  functions  and  some  of  the  nonstandard  ones.  This 
package  has  been  interfaced  with  the  AUGMENT  precompiler  (see  (4]). 

The  multiple  precision  version  of  the  Triplex  Package  is  based  on  Brent's 
multiple  precision  arithmetic  package*  The  multiple  precision  package  is  used 
as  type  EXTENDED,  and  BPA  routines  have  been  written  for  the  primitive  opera- 
tions, using  an  arithmetic  format  which  corresponds  quite  closelv  with  the 
format  of  Brent's  multiple  precision  numbers.  The  basic  architecture  of  the 
Triplex  Package  is,  however,  unchanged. 

In  order  to  use  the  multiple  precision  version,  one  must  have  access  to 
Brent's  package.  This  may  be  obtained  by  writing  directly  to  Brent  at  the 
Computing  Research  Group,  Australian  National  University,  Canberra,  Australia. 

Modification  of  the  single  precision  version  of  the  triplex  package  to 
produce  the  multiple  percision  version  is  entirely  analogous  to  modification 
of  the  single  precision  version  of  the  interval  arithmetic  package  to  produce 
the  multiple  precision  version.  The  latter  process  is  discussed  in  detail  in 
[14],  and  wp  will  not  iterate  the  discussion  here. 

The  basic  multiple  precision  version  of  the  triplex  package  carries  about 
28  decimal  digits,  but  this  is  easily  modified  by  a natural  extension  of  the 
procedure  described  in  [19]. 
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CanqlysiQp: 


We  have  attempted  to  provide  user  documentation  for  the  Triplex  Package, 
together  with  enough  technical  documentation  to  assist  the  user  in  making  mod- 
ifications that  might  be  desirable  in  a production  environment*  The  balance 
of  the  technical  documentation  will,  as  already  noted,  be  found  in  [I7].  Any 
errors  or  inadequacies  in  this  document  should  be  reDorted  to  the  second  au- 
thor* 


Since  both  the  package  and  this  document  are  based  on  the  interval  arith- 
metic package  described  in  [I7],  those  who  contributed  materially  to  that  ef- 
fort are,  by  transitivity,  also  contributors  to  this  package*  They  include 
Dr*  F.  D.  Crary  and  Dr.  S«  T.  Jones,  both  formerly  of  the  Mathematics  Research 
Center;  Mr*  Bill  Boyt  and  Mr.  James  B.  Cheek,  of  the  U.  S.  Army  Corps  of  Engi- 
neers Waterways  Experiment  Station,  Vicksburg,  Mississippi;  Prof.  Bruce  D. 
Shriver  of  the  University  of  Southwest  Louisiana;  Prof.  Ronnie  Ward  of  the 
University  of  Texas  at  Arlington;  Prof.  Richard  Hetherington  of  the  University 
of  Kansas;  and  Prof*  Myron  Ginsberg  of  Southern  Methodist  University. 


DISCLAIMER 

This  program  has  been  tested  and  is  believed  to  be  correct.  However,  in 
any  program  of  this  size,  residual  errors  are  likely*  Neither  the  authors, 
nor  the  Mathematics  Research  Center,  nor  the  University  of  Wisconsin,  nor  the 
United  States  Army,  nor  any  other  party,  conceivable  or  inconceivable,  will 
assume  any  responsibility  for  damages,  whether  direct,  indirect, 
consequential,  or  inconsequential,  arising  from  errors,  omissions,  malfunc- 
tions, irregularities,  inefficiencies,  or  any  other  sins  of  omission  or 
commission  in  this  program  and/or  its  documentation. 
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APPENDIX  1 


STANDARD  FORTRAN  NUMBER  AND 
TRIPLEX  NUMBER  REPRESENTATIONS 


digit 

: : ; 

0!1!2i3!^!5l6|7}8!9 

SIGN 

: : = 

+ i- 

INTEGER 

: : = 

NULL I <SIGN> ! <INTEGERXDIGIT> 

RADIX 

: : = 

• 

FIXEDPOINT 

: : = 

<INTEGERXRADIX>  | <FIXEDPOINTXDIGIT> 

EXPSEP 

• : : = 

E|D 

EXPONENT 

: : = 

<SIGN> <EXPSEP>  | <EXPSEPXSIGN>  I <EXPONENTXDIGIT> 

NUMBER 

: : = 

<INTEGER>  | <FIXEDPOINT>  | <INTEGERXEXPONENT> ! 

<FIXED POINT ><EXPONENT> 

ENDPTSEP 

• • V 

• • • 

• 

• 

COMMA 

• la 

f 

TRIPLEX 

• • 

<NUMBER>|  (<NUMBER>)  1 <NUMBERXENDPTSEPXNUMBER> ! 

(<NUMBERXENDPTSEPXNUMBER>)  | ( <NUMBERXCOMMAXNUMBER> ) 
<NUMBERXENDPTSEPXNUMBERXENDPTSEPXNUMBER>| 
(<NUMBERXENDPTSEPXNUMBERXENDPTSEP><NUMBER>)  ! 
(<NUMBERXCOMMAXNUMBERXCOMMAXNUMBER>) 
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APPENDIX  2 


PACKAGE  MODULES 


This  appendix  contains  information  on  the  modules  of  the  triplex  package* 
The  BPA  portion  of  the  package  is  not  listed  here;  the  interested  reader  is 
instead  referred  to  [I7].  The  following  three  pases  list  the  modules  of  the 
triplex  portion  of  the  package. 

The  first  column  of  the  table  gives  the  function  of  the  module;  the  sec- 
ond column  gives  an  expanded  explanation  of  the  function  or,  in  some  cases, 
the  definition  of  the  function.  The  third  column  gives  the  data  type  of  the 
result  of  the  module,  according  to  the  code  listed  below. 

The  fourth  column  gives  an  example  of  how  the  module  might  be  invoked  in 
a program  which  is  to  be  processed  by  the  AUGMENT  precompiler.  The  reference 
name  of  the  module  (in  AUGMENT)  is  either  the  function  name  or  the  operator 
shown.  Note  that  all  conversion  modules  may  also  be  invoked  implicitly  bv  a 
replacement  operator. 

The  fifth  column  indicates  how  the  module  would  be  invoked  in  the  absence 
of  AUGMENT.  If  the  routine  is  a subroutine  (indicated  by  an  'S'  in  Column  6), 
the  expression  given  must  be  preceded  by  the  word  'CALL';  otherwise,  the  ex- 
pression is  complete.  Note  that  functions  must  be  typed  appropriately  in  the 
calling  program;  Column  6 gives  the  necessary  information  in  these  cases. 

Data  types ; Data  types  are  indicated  by  one-  or  two-letter  abbrevia- 

tions. The  abbreviations  used  are  as  follows: 

B - BPA 

D - DOUBLE  PRECISION 
E - EXTENDED 
H - HOLLERITH  (unpacked) 

I - INTEGER 

IA  - INTEGER  ARRAY 

L - LOGICAL 

PH  - PACKED  HOLLERITH 

R - REAL 

X - TRIPLEX 

Routine  types : S indicates  subroutine;  any  other  letter  indicates  a 
function  of  the  indicated  type  (see  above).  An  asterisk  denotes  a primitive 
routine. 

Variable  names : The  first  letter  (or,  sometimes,  pair  of  letters)  indi- 
cates the  type  of  the  variable.  The  last  letter  is  R for  result,  A,  B,  or  C 

for  argument.  Other  letters  may  sometimes  be  used  for  special  purposes.  Ex- 
planations not  given  here  will  be  found  in  the  text  or,  in  cases  where  the 

routine  is  not  normally  directly  accessible  to  the  user,  in  tne  program  list- 

ing. 
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RESULT  ROUTINE  INVOCATION  ROUTINE 

OPERATION  DEFIN7.T  ION/EXPLANATION  TYPE  VIA  AUGMENT  DIRECT  TYPE 

ARITHMETIC 

Add  Sum  of  two  triplexes  X XA  ♦ XB  TPXADD(XA,  XB,  XR)  S 

Subtract  Difference  of  two  triplexes  X XA  - XB  TPXSUB(XA,  XB,  XR)  S 
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RESULT  ROUTINE  INVOCATION  ROUTINE 

OPERATION  DEFINITION/EXPLANATION  TYPE  VIA  AUGMENT  DIRECT  TYPE 

SPECIAL  TRIPLEX  FUNCTIONS  (Continued) 

Magnitude  Sup(Aba(XA))  B MAG(XA)  TPXMAGfXA,  BR)  S 
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APPENDIX  3 


FAULT  INFORMATION 

FAULT 

NUMBER  MEANING  ACTION 

0 NO  FAULT  0 

1-63  ARITHMETIC  FAULTS  IN  TRIPLEX  NUMBER: 

FAULT  CODE  IS 

1 6#FAULT ( L )+4  *F  AULT ( M) +F  AULT ( R ) 

WHERE  L=LEFT,  M=MAIN,  AND  R=RIGHT, 

AND  FAULT  CAN  HAVE  ANY  OF  THE  FOLLOWING 
VALUES: 

0 - SUCCESSFUL  OPERATION 

1 - OVERFLOW 

2 - INFINITY 

3 - UNDERFLOW 


FOR  EXAMPLE,  A FAULT  VALUE  OF  43  = 

2* 1 6 ♦ 2»4  + 3 WOULD  MEAN  THAT  THE 
LEFT  ENDPOINT  AND  MAIN  VALUE  WERE  BOTH 
(MINUS)  INFINITY,  AND  THE  RIGHT  ENDPOINT 
UNDERFLOWED. 


THE  ACTION  CODES  FOR  THESE  FAULTS  ARE 
SET  AS  FOLLOWS: 

LOGICALLY  IMPOSSIBLE  FAULTS  4 

ANY  FAULT  INVOLVING  INFINITY 

WHICH  IS  LOGICALLY  POSSIBLE  3 

ANY  FAULT  INVOLVING  ONLY  SUCCESSFUL 

OPERATIONS  OR  UNDERFLOWS  0 

64  DIVISION  BY  ZERO  3 

65  ZERO  ••  ZERO  1 

66  SQUARE  ROOT  OF  NEGATIVE  NUMBER  3 

67  LOG  OF  NON POSITIVE  NUMBER  3 

68  UNDERFLOW  DURING  BPA  COMPUTATION  0 

69  OVERFLOW  DURING  BPA  COMPUTATION  3 

70  INTERSECTION  OF  DISJOINT  TRIPLEXES  3 

71  ARC  COS  OR  ARC  SIN  ARGUMENT  OUT  OF  RANGE  3 

72  INVERTED  TRIPLEX  4 

73  ILLEGAL  INPUT  CHARACTER  4 

74  ILLEGAL  INPUT  FORMAT  SPECIFICATION  4 

75  ILLEGAL  OUTPUT  FORMAT  SPECIFICATION  1 

76  INPUT  STRING  TOO  LONG  4 

77  UNSPECIFIED  INPUT  UNIT  4 

78  END  OF  FILE  ON  INPUT  1 

79  UNSPECIFIED  OUTPUT  UNIT  1 

80  CONVERSION  ARRAY  OVERFLOW  4(A) 

81  UNRECOGNIZED  ERROR  4 


(A)  THIS  FAULT  VALUE  SHOULD  NOT  BE  CHANGED,  SINCE  ANY  OTHER 
ACTION  COULD  RESULT  IN  A RECURSIVE  CALL  ON  TPXRAP  FROM 
TPXCXH. 
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TPXRAP  WILL  TAKE  ONE  OF  SEVERAL  POSSIBLE  ACTIONS,  DEPENDING 
ON  THE  VALUE  OF  THE  CELL  IN  THE  MON  ARRAY  CORRESPONDING 
TO  THE  FAULT  WHICH  HAS  OCCURRED,  POSSIBLE  VALUES  AND 
CORRESPONDING  ACTIONS  ARE 

0 RETURN  WITHOUT  TAKING  ACTION 

1 PRINT  MESSAGE 

2 PRINT  MESSAGE  AND  TRACE  CALL  SEQUENCE 

3 AS  (2),  AND  STEP  ERROR  COUNTER 

4 AS  (2),  BUT  THEN  TERMINATE  JOB 

IF  NO  INFORMATION  TO  THE  CONTRARY  IS  SUPPLIED,  TPXRAP 
WILL  ASSUME  THAT  THE  CALLING  ROUTINE  IS  A SUBROUTINE 
WITH  TWO  TRIPLEX  ARGUMENTS  ~ AN  OPERAND  AND  A RESULT, 
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APPENDIX  4 


COMMON  BLOCK  INFORMATION 


This  appendix  contains  information  concerning  the  COMMON  blocks  used  in 
the  Triplex  Package*  This  is  an  abbreviated  version  of  the  corresponding  ap- 
pendix in  [$Y1],  and  the  reader  needing  more  information  should  consult  that 
document.  The  changes  indicated  in  Section  5 will  render  that  appendix  valid 
for  the  Triplex  package. 

Here,  we  list  only  the  COMMON  block  declarations,  together  with  DATA 
statements  for  the  portions  of  the  Triplex  Package  which  are  not  dependent  on 
the  representation  of  BP A,  EXTENDED,  or  TRIPLEX.  Thus,  the  reader  can  see 
what  declarations  are  necessary  and  the  values  assigned  to  the  variables  in 
question.  This  should  enable  one  to  make  any  modifications  that  would  be  de- 
sirable in  any  routine  use  of  the  package. 


1.  COMMON -BLOCK  FOR  DEFINING  ROUNDING  OPTIONS  AND  PASSING 

C SELECTED  ROUNDING  OPTION,  ACCURACY,  AND  RESULTING  FAULTS  TO 

C AND  FROM  THE  BPA  PACKAGE 

COMMON  /BPACOM/  OPTION,  BPAFLT,  ACC,  RDU,  RDL,  RDT , RDN,  RDA 
INTEGER  OPTION,  BPAFLT,  ACC,  RDU,  RDL,  RDT,  RDN,  RDA 
DATA  RDU  /I/,  RDL  /2/,  RDT  /3/,  RDN  /4/,  RDA  /5/ 


2.  COMMON  BLOCK  NEEDED  TO  HOLD  PARAMETERS  FOR  BPA  ROUTINES. 

C THESE  CONSTANTS  ARE  MACHINE  - DEPENDENT. 

COMMON  /BPACON/  ZRO,  ONE,  ONM,  TWO,  PI02,  PI, 


BPA 

DATA  ZRO 


1 

2 

3 

4 

5 

6 

7 


PI, 


BPAMNB,  BPAMXB 
BPAMNB,  BPAMXB 


ZRO,  ONE,  ONM,  TWO,  PI02, 

/<BPA  representation  of  0>/, 

/<BPA  representation  of  1>/, 

/<BPA  representation  of  — 1 > / , 

/<BPA  representation  of  TWO>/, 

/<BPA  representation  of  PI/2,  rounded  ud>/, 

/<BPA  representation  of  PI,  rounded  up>/, 

BPAMNB  /<BPA  representation  of  smallest  positive  BPA  number>/ 
BPAMXB  /<BPA  representation  of  largest  positive  BPA  number>/ 


ONE 

ONM 

2 

PI02 

PI 


3.  COMMON  BLOCK  FOR  CONSTANTS  DEALING  WITH  ACCURACY  OF 

LIBRARY  FUNCTIONS 

COMMON  /BPAACC/  IACC(24),  LNACC,  LGACC,  SNHACC.TNHACC,  EXPMXA, 

1 EXPMNA,  FRACBD , MAXINT 

INTEGER  I ACC 

BPA  LNACC,  LGACC,  SNHACC,  TNHACC,  EXPMXA,  EXPMNA,  FRACBD,  MAXINT 
C VALUES  ARE  GIVEN  FOR  ACCURACY  (IN  BITS)  OF  THE  EXTENDED 

C LIBRARY  FUNCTIONS  SIN,  COS,  TAN,  ASN,  ACS,  ATN,  LN( AWAY  FROM 

C 1),  LN(NEAR  1),  L0G10( AWAY  FROM  1),  LOG10(NEAR  1),  EXP, 

C SNH( AWAY  FROM  0),  SNH(NEAR  0),  COSH,  TANH(AWAY  FROM  0), 

C TANH(NEAR  0),  SORT,  AND  CBRT , RESPECTIVELY. 

DATA  IACC  /<values  for  accuracy  in  digits  of  individual  routines,  as 
1 listed  above/ 

DATA  LNACC  /Cradius  of  nbd  of  1 in  which  LN  is  less  accurate>/, 

1 LGACC  /<radius  of  nbd  of  1 in  which  L0G10  is  less  accurate>/, 
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4 0<t/l  CUI  M 


SNHACC  /<radius  of  nbd  of  0 in  which  SINH  is  less  aceurate>/, 
TNHACC  /<radius  of  nbd  of  0 in  which  TANH  is  less  accurate?/, 
EXPMXA  /< largest  BPA  such  that  EXP(EXPMXA)  does  not  overflow?/, 
EXPMNA  /{smallest  BPA  such  that  EXP(EXPMNA)  does  not  underflow?/, 
FRACBD  /{smallest  BPA  number  with  all  digits  preceding  radix?/, 
MAXINT  /{BPA  representation  of  (largest  integer  •*-  1)?/ 


4.  COMMON  BLOCK  FOR  PASSING  INFORMATION  NEEDED  BY  BPA 

C CONVERSION  ROUTINES 

COMON  /BPACCM/  EXWDTH,  ESDMAX,  EXMAX,  EXMIN, 


1 

IX(80) , ICX(40) , ECX(40), 

EX(80) 

INTEGER 

EXWDTH,  ESDMAX,  EXMAX,  EXMIN, 

1 

IX,  ICX,  ECX,  EX 

DATA 

EXWDTH  121 , ESDMAX  /60/,  EXMAX  /40/, 

EXMIN 

DATA 

IX 

/80, 

2,  1,12,  0,  0,  0,  0,  27/, 

1 

ICX 

/4D, 

65536,  16,  21,  0,  0,  0,  0, 

3/, 

2 

ECX 

/40, 

100000,  5,  21,  0,  0,  0/, 

3 

EX 

/80, 

10,  1,12/ 

5.  COMMON  BLOCK  FOR  COMMUNICATION  WITH  ERROR-HANDLING  ROUTINE 

COMMON  /TPXCOM/  TPXFLT,  ID,  TA,  TB,  TR 
INTEGER  TPXFLT,  ID 
TRIPLEX  TA,  TB,  TR 
EXTENDED  TAE,  TBE,  TRE 

EQUIVALENCE  (TA,  TAE),  (TB,  TBE),  (TR,  TRE) 


6. 


MACHINE  DEPENDENT  TRIPLEX  CONSTANTS  AND  EXTENDED  CONSTANTS 
FOR  REDUCTION  OF  TRIG  FUNCTION  ARGUMENTS  TO  PRINCIPAL  RANGE 
COMMON  /TPXCON/  PI02I,  PII,  TPI02I , TPII,  NINSTH, 

1 EPI,  ETPI,  DELTA,  EONE 

TRIPLEX  PI02I , PII,  TPI02I , TPII,  NIN8TH 

EXTENDED  EPI,  ETPI,  DELTA,  EONE 

REAL  REPI(2),  RETPI(2),  RDELTA(2) , RE0NE(2 ) 

EQUIVALENCE (REPI ( 1 ) , EPI),  (RETPI(I),  ETPI),  (RDELTA(I),  DELTA), 
1 (REONE(I),  EONE) 

DATA  PI02I  /{triplex  representation  of  PI/2?/, 

1 PII  /{triplex  representation  of  PI?/, 

2 TPI02I  /{triplex  representation  of  3PI/2?/, 

3 TPII  /{triplex  representation  of  2PI?/, 

4 NIN8TH  /{triplex  representation  of  9/8?/, 

5 REPI  /{extended  representation  of  PI?/, 

6 RETPI  /{extended  representation  of  2PI?/, 

7 RDELTA  /{precision  of  extended  arithmetic?/, 

8 REONE  /{extended  representation  of  1»0?/ 


7.  CONSTANTS  FOR  TRIPLEX  I/O  ROUTINES 

COMMON  /TPXCCM/  NUNITS,  UNITBL(8,  5),  NOCHR,  ICHR(IO),  OCHR(4), 

1 LRBFOR , RBF0RMO32),  LRBFRE,  RBFREE(  132),  IRBFRE,  OLDUNT, 

2 LWBFOR,  WBFORM(132),  SOFMT(4),  IWMIN,  NOCCC,  CCC(2), 

3 KFMT(4 ) , HSTRH32) 

INTEGER  NUNITS,  UNITBL,  NOCHR,  ICHR,  OCHR,  LRBFOR,  RBFORM,  LRBFRE, 

1 RBFREE,  IRBFRE,  OLDUNT,  LWBFOR,  WBFORM,  SOFMT , IWMIN, 

2 NOCCC,  CCC,  KFKT,  HSTR 
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DATA  NUNITS  /4/,  ((UNITBL(I,  J ) , J = 1 , 5) , I = 1 , 4) 

1 / 5,  80,  80,  0,  1, 

2 6,  132,  132,  1,  2, 

3 1,  80,  80,  0,  2, 

4 0,  80,  80,  0,  1/ 

DATA  NOCHR  /10/ 

DATA  ICHR  / 1 H ( , 1H:,  1H),  1H  , 1H,,  1H#,  1H=,  1H$,  1H0,  1H  / 

DATA  OCHR  /1H(,  1H,,  1H),  1H  / 

DATA  IBUFRE  / 1 33/ , OLDUNT  /-I/,  LRBFOR  /1 32/,  LRBFRE  /1 32/ , 

1 LWBFOR  /1 32/,  SOFMT  /2,  1,  55,  1H  /,  IWMIN  /22/,  NOCCC  /2/, 

2 CCC  /1H  , 1H0/ 

COMMON  BLOCK  TO  CONTROL  ECHOING  OF  INPUT  RECORDS  AND  FIELDS 
COMMON  /TPXCRD/  ECHOS,  ECHOD 
LOGICAL  ECHOS,  ECHOD 
DATA  ECHOS  /.FALSE./,  ECHOD  /.FALSE./ 

COMMON  BLOCK  DECLARATIONS  FOR  TPXRAP 

THE  PURPOSE  OF  THIS  COMMON  BLOCK  IS  TO  GIVE  THE  USER 
ACCESS  TO  THE  PARAMETERS  USED  BY  TPXRAP  WITHOUT  REQUIRING 
RECOMPILATION  OF  TPXRAP. 

COMMON  /TPXRCM/  MON(88),  NAMES(40),  TYPA(40),  TYPB(40),  TYPR(40) 
INTEGER  MON,  NAMES,  TYPA,  TYPB , TYPR 


DATA 

MON 

/0,  4,  3, 

0,  4,  4,  3,  4,  4,  4, 

3,  4,  0,  4,  3,  0, 

1 

*, 

*».  *»,  *».  3,  4,  4,  4, 

3,  4,  4,  4,  4,  4, 

2 

3,  *»,  3, 

3,  3,  3,  3,  3,  3,  3, 

3,  3,  3,  4,  3,  3, 

3 

o,  4,  3, 

0,  4,  4,  3,  4,  4,  4, 

3,  4,  0,  4,  3,  0, 

4 

3,  1,  3, 

3,  0,  3,  3,  3,  4,  4, 

4,  1,  4,  4,  1,  1,  4, 

4/ 

DATA 

NAMES  /3HADD, 

3HSUB,  3HMUL,  3HDIV, 

3HSCT,  3HXXI,  3HSIN, 

1 

3HCOS, 

3HTAN,  3HASN,  3HACS, 

3HATN,  3HEXP,  3HLN  , 

2 

3HLOG, 

3HSNH,  3HCSH,  3HTNH, 

3HSQT,  3HCBT,  3HMDB, 

3 

3HHLB, 

3HLGB,  3HCXI,  3HBND, 

3HCEX,  3HRD  , 3HRDF, 

4 

3HASG, 

3HCHX,  3HWR  , 3HCXH, 

3HCRX,  3HDST,  3HCPS/ 

DATA 

TYPA 

/6,  6,  6, 

6,  6,  6,  6,  6,  6,  6, 

6,  6,  6,  6,  6,  6,  6, 

6, 

1 

6 , 6 , 6 , 

6,  6,  6,  5,  5,  0,  0, 

0,  0,  0,  0,  2,  6,  0/ 

DATA 

TYPB 

/6 , 6,  6, 

6,  6,  1,  0,  0,  0,  0, 

0,  0,  0,  0,  0,  0,  0, 

o, 

1 

0 , 0 , 0 , 

0,  0,  0,  1,  0,  0,  0, 

0,  1,  0,  1,  0,  6,  0/ 

DATA 

TYPR 

/6,  6,  6, 

6,  6,  6,  6,  6,  6,  6, 

8,  6,  8,  6,  6,  6,  8, 

6, 

1 

6,  6,  4, 

^ 6,  0,  0, 

0,  6,  0,  6,  6,  4,  6/ 

i.wiiwwyw. 
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TRIPLEX  DESCRIPTION  DECKS  FOR  AUGMENT 


DSTOOIOO 
DST001 10 


SINGLE  PRECISION  VERSION  OF  TRIPLEX  PACKAGE 

•DESCRIBE  EXTENDED  DST00010 

DECURE  DOUBLE  PRECISION,  KIND  FUNCTION,  PREFIX  EXT  DST00020 

COMMENT  DUMMY  DESCRIPTION  SO  CONVERT  CAN  BE  USED  DST00030 

FUNCTION  INT  (=DINT (STD) , (DOUBLE  PRECISION),  DOUBLE  PRECISION)  DST00040 

•DESCRIBE  BPA  DST00050 

DECURE  REAL,  KIND  SAFE  SUBROUTINE,  PREFIX  BPA  DST00060 

OPERATOR  + (,  NULL  UNARY,  PRV,  $)  DST00070 

COMMENT  THE  FOLLOWING  DESCRIPTOR  MUST  BE  INSERTED  IN  THE  DECK  AT  THIS  DST00080 
POINT  IF  THE  REAL  UNARY  MINUS  OF  FORTRAN  DOES  NOT  SUFFICE  FOR  THE  BPA  DST00090 
UNARY  MINUS  — DSTOOIOO 

OPERATOR  - (NEG,  UNARY,  PRV,  $)  DST00110 

OPERATOR  + (ADD,  BINARY  1,  PRV,  $,  COMM),  • (MUL),  - (SUB, , , , NONCOMM) , DST00120 
/ (DIV)  DST00130 

OPERATOR  ••  (XBI,  BINARY  3,  PRV,  $,  INTEGER,  $)  DST00140 

OPERATOR  .EQ.  (EQ,  BINARY  2,  PRV,  $,  LOGICAL),  .NE.  (NE),  .LT.  (LT),  DST00150 

.LE.  (LE),  .GT.  (GT),  .GE.  (GE)  DST00160 

FUNCTION  ABS  (ABS,  ($),  $),  SIN  (SIN),  COS  (COS),  TAN  (TAN),  ASIN  (ASN) ,DST00170 
ACOS  (ACS),  ATAN  (ATN) , L0G10  (LOG),  LOG  (LN),  LN  ( ) , DST00180 

EXP  (EXP),  SINH  (SNH) , COSH  (CSH),  TANH  (TNH) , SORT  (SOT),  DST00190 
CBRT  (CBT)  DST00200 

FUNCTION  INT  (INT,  ($),  $)  DST00210 

FUNCTION  MAX  (MAX,  ($,  $),  $),  MIN  (MIN)  DST00220 

FUNCTION  SGN  (SGN , ($),  INTEGER)  DST00230 

COMMENT  THE  FOLLOWING  DESCRIPTOR  MUST  BE  INSERTED  IN  THE  DECK  AT  THIS  DST00240 
POINT  IF  THE  REAL  REPLACEMENT  OPERATOR  IN  FORTRAN  WILL  NOT  SERVE  AS  THEDST00250 
BPA  REPLACEMENT  OPERATOR  DST00260 

SERVICE  COPY(STR)  DST00270 

CONVERSION  CTB  (CEB,  DOUBLE  PRECISION,  $,  DOWNWARD),  DST00280 

CTB  (CRB,  REAL,  $,  UPWARD),  DST00290 

CTB  (CIB,  INTEGER,  $,  UPWARD),  DST00300 

CTE  (CBE,  $,  DOUBLE  PRECISION,  UPWARD),  DST00310 

CTR  (CBR,  $,  REAL,  DOWNWARD),  DST00320 

CTI  ( CBI , $,  INTEGER,  DOWNWARD)  DST00330 

COMMENT  THE  FOLLOWING  TWO  CARDS  ARE  USED  TO  FORCE  AUGMENT  TO  USE  THE  DST00340 
REAL  UNARY  MINUS  OPERATOR  OF  FORTRAN  AS  THE  BPA  UNARY  MINUS  OPERATOR  — DST00350 
•ENVIRONMENT  DST00360 

OPERATOR  - (UNARY,  PRV,  BPA)  DST00370 

•DESCRIBE  TRIPLEX  DST00380 

COMMENT  IN  THIS  DESCRIPTION  DECK,  'DOUBLE  PRECISION*  IS  USED  AS  A DST00390 

SYNONYM  FOR  'EXTENDED'.  IF  EXTENDED  IS  A TYPE  OTHER  THAN  DOUBLE  PRE-  DST00400 
CISION,  THE  APPROPRIATE  CHANGES  WILL  NEED  TO  BE  MADE  IN  THIS  DECK.  DST00410 
DECURE  REAL(3),  KIND  SAFE  SUBROUTINE,  PREFIX  TPX  DST00420 

OPERATOR  ♦ (.NULL  UNARY,  PRV,  $)  DST00430 

OPERATOR  - (NEG,  UNARY,  PRV,  $)  DST00440 

OPERATOR  ♦ (ADD,  BINARY1 , PRV,  t,  COMM),  • (MUL),  - (SUB, ,,, NONCOMM) , DST004S0 

/ (DIV)  DST00460 


P 

P 


OPERATOR  »»  (XXI,  BINARY  3,  PRV,  $,  INTEGER,  $),  »•  (XXB, , , ,BPA) , DST00470 

*•  (XXE, ,, .DOUBLE  PRECISION),  ••  (XXX,,,, $)  DST00U80 

OPERATOR  .VEO.  (VEQ,  BINARY  2,  .EQ.  , $,  LOGICAL,  COMM),  .VNE.  (VNE),  DST00490 
.SEQ.  (SEQ),  .SNE.  (SNE)  DST00500 

OPERATOR  .VLT.  (VLT,  BINARY  2,  .LT.  , $,  LOGICAL),  .VLE.  (VLE) , DST00510 

.VGT.  (VGT ) , .VGE.  (VGE) , .SLT.  ( SLT ) , .SLE.  (SLE),  DST00520 

.SGT.  (SGT),  .SGE.  (SGE ) DST00530 

OPERATOR  .UNION.  (UNN,  BINARY  1,  .OR.,  $),  .INTSCT.  (SCT,,  .AND.)  DST00540 

OPERATOR  .SUBSET.  (SBS,  BINARY  2,  .EQ.,  $,  LOGICAL),  .SPRSET.  (SPS)  DST00550 

OPERATOR  .E.  (ELE,  BINARY  3,  .EQ.,  BPA,  $,  LOGICAL)  DST00560 

FUNCTION  ABS  (ABS,  ($),  $),  SIN(SIN),  COS  (COS),  TAN  (TAN),  ASIN  (ASN) , DST00570 
ACOS  (ACS),  AT AN  ( ATN) , LOG 10  (LOG),  LN  (LN),  EXP  (EXP),  DST00580 

SINH  (SNH) , COSH  (CSH),  TANH  (TNH) , SQRT  (SOT),  CBRT(CBT),  DST00590 

LOG  (LN)  DST00600 

FUNCTION  ATAN2  ( AT2,  ($,  $),  $)  DST00610 

FUNCTION  OK  (OK,  ($),  LOGICAL),  BAD  (BAD)  DST00620 

FUNCTION  COMPOS  (CPS,  (BPA,  BPA,  BPA),  $)  DST00630 

FUNCTION  INT  (INT,  ($),  $)  DST00640 

FUNCTION  DIST  (-OST,  ($,  $),  BPA)  DST00650 

FUNCTION  MDPT  (MDB,  ($),  BPA),  HLGTH  (HLB) , LENGTH  (LGB) , LGTH  (LGB) , DST00660 

MAG  (MAG),  MIG  (MIG),  PIVL  (PVL),  PIVU  (PVU),  SIZE  (SIZ)  DST00670 

FUNCTION  SGN  (SGN,  ($),  INTEGER)  DST00680 

FUNCTION  BOUND  (BND,  (DOUBLE  PRECISION,  INTEGER) , $)  DST00690 

FIELD  SUP  (SUP,  SPL,  ($),  BPA),  INF  (INF,  INL) , MAIN  (MNR , MNL)  DST00700 

SERVICE  COPY  (STR)  DST00710 

CONVERSION  CTX  (CIX,  INTEGER,  $,  UPWARD),  DST00720 

CTX  (CBX,  BPA,  $,  UPWARD),  DST00730 

CTX  (CRX,  REAL,  $,  UPWARD),  DSTOO^tO 

CTX  ( CEX,  DOUBLE  PRECISION,  $,  UPWARD),  DST00750 

CTX  (ASG , HOLLERITH,  $,  UPWARD),  DST00760 

CTI  (CXI,  $,  INTEGER,  DOWNWARD),  DST00770 

CTB  (CXB,  $,  BPA,  DOWNWARD),  DST00780 

CTR  (CXR,  $,  REAL,  DOWNWARD),  DST00790 

CTE  (CXE,  $,  DOUBLE  PRECISION,  DOWNWARD)  DST00800 
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APPENDIX  6 


RUNSTREAM  FOR  USE  OF  THE  TRIPLEX  PACKAGE  WITH  AUGMENT 

We  first  give  an  outline  of  the  AUGMENT  runstream  which  should  be  appli- 
cable to  and  AUGMENT  installation. 

1.  [ATTACH  APPROPRIATE  FILES  TO  OBTAIN  AUGMENT  AND 

THE  DESIRED  NONSTANDARD  PACKAGES] 

2.  [INVOKE  AUGMENT] 

PRINT  SUPPRESS  (if  listing  of  the  description  decks 
is  not  desired) 

3.  [SUPPLY  APPROPRIATE  DESCRIPTION  DECKS  TO  AUGMENT] 

•BEGIN 

<any  special  cards> 

•DISABLE  WARNINGS  (if  AUGMENT  warning  messages 

are  not  desired) 

4.  [SOURCE  CODE  TO  BE  PROCESSED  BY  AUGMENT] 

•END  - 

5.  [COMPILE  PREPROCESSED  PROGRAM,  WHICH  HAS  BEEN  WRITTEN 

ON  LOGICAL  UNIT  20] 

6.  <any  other  processing,  such  as  preparation  of  files,  compi- 

lation of  other  modules,  etc.> 

7.  [INVOKE  LINKAGE  EDITOR] 

<special  instructions  to  linkage  editor,  aside  from 
the  rather  obvious  requirement  that  the  nonstandard 
package  must  be  loaded  with  the  program> 

8.  [EXECUTE  PROGRAM] 

<data> 

9.  [END  OF  JOB] 

The  following  runstream  for  the  single  precision  version  of  the  Tridex 
Package  on  the  UNIVAC  1110  installation  at  the  University  of  Wisconsin  - Madi- 
son shows  how  the  above  runstream  is  implemented  in  a specific  situation: 

§ASG , A MRC*LIB. 
euSE  A. ,MRC«LIB. 

§ASG , A MRC»TRIPLEX. 

#USE  P.,MRC»TRIPLEX. 

§XQT  A. AUGMENT 

PRINT  SUPPRESS  (if  listing  of  the  description  decks 
is  not  desired) 

0ADD  A .DESCRIPTION /FORTR AN- V 
0ADD  P. DESCRIPT ION /TRIPLEX 
•BEGIN 

•CONVERT  EXTENDED  - DOUBLE  PRECISION 
•DISABLE  WARNINGS  (if  AUGMENT  warning  messages 

are  not  desired) 

:F0R, <options>  <filename ,etc.>  (note  rather  than 

n@"  in  Column  1) 

<program> 

: FOR, etc. 

<subprogram> 
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r 


etc* 

•END 
0ADD  20. 

<any  other  processing,  such  as  0PREP  of  files,  compi- 
lation of  other  modules,  etc*> 

£MAP,I  <absolute  element> 

IN  TPF$. 

IN  A.BPACOM 
IN  P.TPXCCOM 
LIB  A. 

LIB  P. 

@XQT  <absolute  element > 

<data> 

@FIN 
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